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We solve analytically the N x N square lattice Hubbard model for even at half 
filling and weak coupling by a new approach. The exact ground state wave function 
provides an intriguing and appealing picture of the antiferromagnetic order. Like at 
strong coupling, the ground state has total momentum Ktot ~ (0, 0) and transforms 
as an s wave for even N/2 and as a d wave otherwise. 



The 2D Hubbard model is one of the simplest descriptions of itinerant interacting electrons on a 
lattice. The Lieb theorem states that at half filling the ground state for a bipartite lattice is unique 

and has spin " where (|52|) is the number of sites in the Si {S2) sublattice; here and in 

the following, \S\ will be the number of elements in the set S. For a, N x N (with even N) square 
lattice, the ground state is a singlet. In Ref. M it was shown that in the strong coupling limit the 
ground state has total momentum Ktot ~ (0, Oj and has s-wave or d-wave symmetry for even or odd 
N/2, respectively. 

In this paper we build the exact ground state of the Hubbard model at half filling in the weak 
coupling limit. We find the same quantum numbers as predicted by Refs. (l) j2|, and this result 
strongly supports the conjecture that no phase transition takes place for finite values of the Coulomb 
interaction parameter U. 

Let us consider the Hubbard model with hamiltonian 

(T {r,7-') r 

on a square lattice of N x N sites with periodic boundary conditions and even N. Here a =t,J, 
is the spin and r, r' the spatial degrees of freedom of the creation and annihilation operators 
and c respectively. We represent sites by r = {ix,iy) with ix,iy = 0,... ,7V — 1. The sum on 
{r,r') is over the pairs of nearest neighbor sites and fira is the number operator on the site r of 
spin a. The point symmetry is C^v, the Group of a square We Fourier expand the fermion 
operators: Cko = ;^ e*'"'cro-, where k = {kx,ky) = ^{ix,iy); then, = Sfe e(^)4<TCfeo- where 
e(fc) = 2t[cos kx+cos ky]. The starting point is the following property of the number operator fir = c\.Cr 
(for the moment we omit the spin index). 

Theorem: Let S he an arbitrary set of plane-wave eigenstates {|fci}} of Hq and {nj.)ij = {ki\fi.r\kj) = 
^Qi-(K-kj)r rnatrix of in S. This matrix has eigenvalues Ai = and A2 = . . . = A|5| — . 

Note that |5| < iV^; if |iS| — iV^ the set is complete, like the set of all orbitals, and the theorem 
is trivial (a particle on site r is the eigenstate with eigenvalue 1). Otherwise, if \S\ < N^, the 
theorem is an immediate consequence of the fact that 



det|(n.).,-A%| = (-A)l'5|-i(M-A), Vr. (2) 



Let Shf denote the set (or shell) of the k wave vectors such that e(fc) = 0. These k vectors lie on 
the square having vertices (±7r, 0) and (0, ±7r); one can show that the number of solutions, that is, the 
dimension of the set Shf, is \Shf \ = 2N — 2. At half filling {N'^ holes) for U = 0, Shf is half occupied. 

/ 2N - 2 \ ^ 

In the Sz = sector there are ( ^ ^ j degenerate unperturbed ground state configurations for 

iV^ holes, where all \k) orbitals such that e(fc) < are occupied. The first-order splitting of the 
degeneracy is obtained by diagonalizing the W matrix over the unperturbed basis; like in elementary 
atomic physics, the filled shells just produce a constant shift of all the eigenvalues and for the rest 
may be ignored in first-order calculations. In other terms, we consider the truncated Hilbert space Ti. 
spanned by the states of — I holes of each spin in Shf, and we want the exact ground state(s) of 
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W in 7i; by construction H is in the kernel of Hq, so the ground state of W is the ground state of 
H as well. We call W — state any vector in H which also belongs to the kernel of W. Since the 
lowest eigenvalue of W is zero, it is evident that any W — state is a ground state of H. We want 
to calculate the unique ground state of the Hubbard Hamiltonian for U — at half filling which is 
a, W — singlet with 2N ~ 2 holes in Shf (filled shells are understood). 

To diagonalize the local operator W in closed form we need to set up a local basis set of one-body 
states. If Shf were the complete set of plane- wave states |fc), the new basis would be trivially obtained 
by a Fourier transformation, but this is not the case. The question is: how can we define for each site 
r the local counterparts of k states using only those that belong to a degenerate level? The answer is: 
build a set {\(pa'')} of orbitals such that the number operator fir and the Dirac characters of the point 
symmetry Group are diagonal. Using such a basis set for the half-filled shell the unique properties 
of the antiferromagnetic ground state become simple and transparent. The eigenvectors l^Pa"^) of nr=o 

\ (r)\ . 

and those \(pa ) of other sites r are connected by translation and also by a unitary transformation, or 
change of basis set. Picking r — ei, I — x means ei — (1,0) or transfer by one step towards the right 
and I = y means ei = (0, 1) or transfer by one step upwards. The unitary transformation reads: 

2N-2 2N~2 
13=1 13=1 

The transfer matrix Ti knows all the translational and point symmetry of the system, and will turn 
out to be very special. 

For large N, to find {ji^si'^'')} it is convenient to separate the fc's oi Shf in irreducible representations 
(irreps) of the space Group [|| dv ® T; here T is the Abelian Group of the translations and (8) 
means the semidirect product. Choosing an arbitrary k € Shf with kx > ky > 0, the set of vectors 
Rik, where Ri £ C4„, is a (translationally invariant) basis for an irrep of G; the accidental degeneracy 
of several irreps is due to the presence of extra symmetry, i.e. G is a subgroup of the Optimal Group 
defined in 0. The high symmetry vectors (0,7r) and (7r,0) always trasform among themselves 
and are the basis of the only two-dimensional irrep of G, which exists for any N. If N/2 is even, one 
also finds the high symmetry wavevectors k = (±7r/2, ±7r/2) which mix among themselves and yield 
a four-dimensional irrep. In general, the vectors Rik are all different, so all the other irreps of G have 
dimension 8, the number of operations of the point Group 6*4^. 

Next, we show how to build our local basis set and derive W = states for each kind of irreps of 
G. For illustration, we will consider the case = 4; then Shf contains the bases of two irreps of G, 
of dimensions 2 and 4. The one with basis fc^ = (tt, 0), fcs = (0, vr) breaks into Ai Bi in C^v ■ 

The eigenstates of {nr=o)ij = {ki\nr=o\kj) , with i,j = A,B , are IV'^^) = -^(I^a) + l^s)) with 
Ai = 1/8 and \ipBi) — ^f^l^^^ ~ l^s)) with A2 = 0. Since under translation by a lattice step T/ along 

the I = x,y direction |fc) e^^'\k), using Equation (||) one finds that \tp'J^^) ^ (~1)^' Wbi)^ with 
6'^ — 1, 9y — 0; so ) has vanishing amplitude on a sublattice and |'0^^ ) on the other. The two-body 
state |?/'^^}(j|'(/'Bj)-cr has occupation for spin a but not for spin —a on the site r = 0; under a lattice 
step translation it flips the spin and picks up a (-1) phase factor: |V'^i)ct|V'bi)-<t ^ |'0Bj)(t|V'ai)-<t; 
therefore it has double occupation nowhere and is a = state (W = pair [|j ||]). 

The 4-dimensional irrep with basis fci = (tt/2,7t/2), k2 = (— 7r/2, 7r/2), k^ ~ (7r/2,— 7r/2) ^4 = 
(-7r/2, -7r/2) breaks into Ai ® B2 ® E in 6*4^; letting /= 1,2,3,4 for the irreps Ai, B2, Ex, Ey 
respectively, we can write down all the eigenvectors of {nr=o)ij = {ki\nr=o\kj) , with i, j — 1, . . . ,4, as 
= J2i=i 0'ji\ki), where O' is the following 4x4 orthogonal matrix 



-4 



1111 

1-1-1 1 
1-1 1-1 
-1-1 1 1 



(4) 



The state with non-vanishing eigenvalue is again the Ai eigenstate. After a little bit of algebra 
we have shown Q that under T; the subspace of Ai and B2 symmetry is exchanged with the one 
of Ex and Ey symmetry. Thus we can build a 4-body eigenstate of W with vanishing eigenvalue: 
WAii^Bi) ''\'^'e V'b )-<t- As before under a lattice step translation this state does not change its spatial 
distribution but a — > —a without any phase factor: IV-'AiV-'bj)^!'*/'^; V'b )-cr ^ IV'b 4''e )o-|V-'^i V'B2)-o"• 
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Now we use these results to diagonalize nr=o on the whole set Shf (we could have done that 
directly by diagonalizing 6x6 matrices but we wanted to show the general method). The eigenstate 
of 71^=0 with nonvanishing eigenvalue always belongs to Ai. The matrix has eigenvalues 3/8 and 
(5 times) 0, as predicted by eq.(||). For r = the eigenvector of occupation 3/8 is {(p^i'') = = 

"TslV'Ai) + \I\Wai)- The other Ai eigenstate of rtr=o has eigenvalue and reads: |02"^) = \<l>2\i) — 



The other eigenvectors , whose symmetry differs from Ai, are 

) = We,)- 



4°^) - 10 



a(o) 



) = and 



(0) 



matrices Ti of Equation (||) such that | ^) = X]j 



6^°))Ti,,,are: 



One finds | 



that the transfer 
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(5) 



The reason why this choice of the basis set is clever is now apparent. The local basis at any site r 

splits into the subsets Sa = {l^i'^li)' \'^2'ai)' I^b])}' ^^'^ = {I'/'sj), I 0b]) 7 I'^e])}' ^ ^^^^^ by a lattice 
step sends members of Sa into linear combinations of the members of Sb, and conversely. 
Consider the 6-body eigenstate of i/o 



(6) 



In this state there is partial occupation of site r = with spin cr, but no double occupation. It turns 
out that a shift by a lattice step produces the transformation 



\^AF)a 



(7) 



that is, a lattice step is equivalent to a spin flip, a feature that we call antiferromagnetic property. 
Since the spin-flipped state is also free of double occupation, \^AF)a is a = eigenstate. A ground 
state which is a single determinant is a quite unusual property for an interacting model like this. 

Note that \4''iAi'f''2Ai) equivalent to IV'aiV-'ai)' because this is just a unitary transformation of 
the Ai wave functions; so |$af)o- can also be written in terms of the old local orbitals (without any 
mix of the local states of different irreps of G): 



(8) 



This form of the ground state lends itself to be generalised (see below). For N > i, k vectors arise 
that do not possess any special symmetry, the vectors Rik are all different for all Ri G C^y, and 
we get eight-dimensional irreps of G. Recalling that \Shf\ — 2N — 2, one finds that Shf contains 
Ne = ^(-y ~ 2) irreps of dimension 8, one of dimension 4 and one of dimension 2 if N/2 is even and 
No — |(y ~ 1) irreps of dimension 8 and one of dimension 2 if N/2 is odd. 

To extend the theory to general N, we note that these k vectors, since Rik are all different, are a 
basis of the regular representation of C^y . Thus, by the Burnside theorem, each of them breaks into 
Ai ® A2 ® Bi ® B2 (B E ® E; diagonalizing 71^=0 and the point Group characters on the basis of the 
m-th eight-dimensional irrep of G one gets one-body states where / stands for the C^y irrep 

label, I=Ai, A2, Bi, B2, E^, Ey, E'^, E'^; here we denote by E' the second occourrence of the irrep 
E. The ground state wave function in Ti. for the half filled case is a generalized version of Equation 
(0). For even N/2, we have 



1$^^). ^ i( n <c'CVL7)^:..^k^:;.).i( n vA7CVgVLt)^kV'^.^B.)-., o) 



m— 1 
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with (7 =t, i- For odd N/2, \^AF)a is similar but the maximum m is No and the jf/;') states do 
not occour. \^AF)cr is a = state, transforms into —\^AF)-a for each lattice step translation 
and manifestly shows an antiferromagnetic order (antiferromagnetic property). Since is a positive 
semidefinite operator \^AF)a is actually a ground state. In the basis of local states these are the only 
two determinantal states (cr =1, |) with the above properties. 

A few further remarks about |$a_f)o- are in order. 1) Introducing the projection operator Ps 
on the spin S subspace, one finds that Ps\<i>AF)a = \'^af)<j 7^ ,¥5 = 0,... , TV - 1. Then, 
a{^AF\W\^AF)a = J2s=i <y (^afI^I^af) a — 0, and this implies that there is at least one ground 
state oi W in Ti, for each S. The actual ground state of H at weak coupling is the singlet |<I>^^)o.. 
2) The existence of this singlet W — ground state is a direct consequence of the Lieb theorem Q. 
Indeed the maximum spin state \^AF^)a is trivially in the kernel of W; since the ground state must 
be a singlet it should be an eigenvector of W with vanishing eigenvalue. 3) The above results and 
Lieb's theorem imply that higher order effects split the ground state multiplet of H and the singlet 
is lowest. 4) The Lieb theorem makes no assumptions concerning the lattice structure; adding the 
ingredient of the G symmetry we are able to explicitly display the wave function at weak coupling. 

Using the expHcit form of Ps=o one finds that Ps=o\^AF)a = —Ps=o\^af)-<7- This identity allows 
us to study how the singlet component transforms under translations, reflections and rotations. We 
have found Q that the total momentum is Ktot — (0,0). To make contact with Ref. |^ we have also 
determined how it transforms under the C^y operations with respect to the center of an arbitrary 
placquette. It turns out that it is even under reflections and transforms as an s wave if N/2 is even 
and as a d wave if N/2 is odd. The present approach lends itself to obtain exact results for other 
fillings as well, but this will be shown elsewhere 
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